Soliton— antisoliton pair production in particle collisions 
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We propose general semiclassical method for computing the probability of soliton-antisoliton pair 
production in particle collisions. The method is illustrated by explicit numerical calculations in 
(1 + l)-dimensional scalar field model. We find that the probability of the process is suppressed by 
an exponentially small factor which is almost constant at high energies. 
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Ever since the discovery of topological solitons, a ques- 
tion remains [H of whether soliton-antisoliton (SA) 
pair can be produced at sizable probability in collision 
of two quantum particles. This process, which involves a 
transition from perturbative two-particle state to a non- 
perturbative state containing SA pair, eludes treatment 
by any of the standard methods. A general expecta- 
tion Q is that the probability of the process is ex- 
ponentially suppressed in weak coupling regime. 
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where E is the total energy, ^ <C 1 is the coupling con- 
stant. Indeed, crudely speaking, one can think of solitons 
as bound states of Ns ^ 1/ particles [l| . Then the sup- 
pression ([1]) is due to multiparticle production [s^. 

In this Letter we propose general semiclassical method 
for computing the leading suppression exponent F{E) of 
the inclusive process "two particles SA pair + par- 
ticles." As a by-product, we calculate the exponent 
Fn{E) of the same process with N initial particles. In 
our method the problem is deformed by introducing a 
small parameter Sp which turns the process of SA pair 
production into a well-known tunneling process. To the 
best of our knowledge, no method of this kind has ever 
been proposed before. 

For definiteness we consider (1 + l)-dimensional scalar 
field theory with action Q 
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This model possesses topological solitons if the scalar po- 
tential V{(j)) has a pair of degenerate minima (j)- and 
see the inset in Fig. [H solid line. Soliton and antisoliton 
solutions interpolate between the minima; their profiles 
are shown in Fig. [2^. 

An obstacle to the semiclassical description of SA pair 
production is related to the fact that soliton and anti- 
soliton attract each other and annihilate classically into 
NsA ^ ^/g^ particles. Thus, there is no potential bar- 
rier separating SA pair from the particle sector and the 
process under study cannot be treated as potential tun- 
neling. 

We get around this obstacle by introducing the po- 
tential barrier between SA pair and perturbative states. 
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FIG. 1. Solutions in (E, N) plane at Sp — 0.4. Numbers near 
the lines with empty and filled points are the values of T and 
^, respectively. Inset: Potential density V^cf)). 



Namely, we add negative energy density {—5p) to the 
vacuum see dashed line in the inset in Fig. [H This 
turns (j)- and 0+ into false and true vacua, respectively; 
the process of SA pair production is now interpreted as 
false vacuum decay [HI induced by particle collisions. The 
latter is a well-studied tunneling process [gI, Q • In the 
end of calculation we will take the limit 5p ^ {). 

The height of the potential barrier between the false 
and true vacua is given by the energy Ech of the crit- 
ical bubble [H — unstable static solution "sitting" on 
top of the barrier. The pressure 5p inside this bubble 
is balanced by the soliton-antisoliton attraction. Let us 
estimate the critical bubble size dch at small 5p. The 
attractive force Fatt between the soliton and antisoliton 
is proportional to the Yukawa exponent exp(— m+(ic6), 
where m\ = V''(0+). Setting Fatt = ^P^ one finds 
deb ^ — log(const • Sp)/m-^. We see that at (5p ^ the 
critical bubble turns into a widely separated SA pair and 
Ecb 271^5, where Ms is the soliton mass. 

Another difficulty is met in the case of = 2 par- 
ticles in the initial state, since states with few quanta 
cannot be described semiclassically. We solve thisprob- 
lem by Rubakov-Son-Tinyakov (RST) method [ll, ^. 
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FIG. 2. (a) Soliton-antisoliton configuration, (b) Contour 
in complex time, (c) — (e) Semiclassical solutions (j)s{t,x) at 
Sp = 0.4; respective values of {E, N) are shown by squares 
with letters in Fig.[Tl In particular, (c) E < Ecb, (d) E ^ Ecb, 
(e) E > Ecb- Color represents arg(/)s. 



Semiclassical method for the calculation of the multi- 
particle probability ^ dit N ^ 1 and ^ <C 1 has been 
developed in Ref. M. This method is based on the saddle- 
point evaluation [13| of the path integral for the proba- 
bility. Below we list the boundary conditions for the 
complex saddle-point configuration 6s{t^x) G C and give 
expression for Fn{E); see Refs. [sl, 0] for derivation. 

Configuration (j)s(t,x) satisfies classical field equation 
SS/6(j) = along the contour in complex time shown in 
Fig. [IJd, where the Euclidean part corresponds to tun- 
neling. In the asymptotic past (j)sit^x) is a collection of 
linear waves above the false vacuum (/>_, 
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as t ^ — oo , (5) 
2 - V'\(t)-). The first boundary 



condition is relation between the amplitudes [8] , 

ak = e-^^^^-%k , (6a) 



where T and are real parameters. In the asymptotic fu- 
ture (j)s{t^x) contains a bubble of true vacuum. The sec- 
ond boundary condition is asymptotic real-valuedness: 



The method is based on the assumption that as long as 
the number of colliding quanta is semiclassically small, 
<C 1/^^, the leading suppression exponent Fn{E) is 
universal, i.e. does not depend on N. In detail, consider 
the A^-particle inclusive probability. 
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where S is S-matrix. The sums in Eq. (|3]) run over all ini- 
tial states of energy E and multiplicity N and final states 
containing SA pair and arbitrary number of quanta. RST 
conjecture states that the multiparticle suppression ex- 
ponent Fn{E) in Eq. (j3j) coincides with the two-particle 
exponent F{E) for N <C l/g^- On the other hand, at 
A^ ^ 1 the initial states in Eq. dU can be described 
semiclassically. Thus, F{E) is computed by evaluating 
semiclassically Fn{E) and taking the small- A" limit 
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It is worth noting that the conjecture (j4]) has been sup- 
ported by field theory calculations [lo| and proved in the 
context of multidimensional quantum mechanics |llLll2|. 

After modifying the original problem we arrived at 
false vacuum decay induced by collision of A" ^ 1 parti- 
cles. The latter process can be described semiclassically, 
as outlined below. The price to pay is the limits 6p ^ 
and g^N 0, which should be taken in the end of calcu- 
lation. In particular, we are going to show that the limit 
Sp ^ of the suppression exponent exists. 
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Equations (|6]) are sufficient to specify complex solution 
(j)s{t^x) for given values of T and in Eq. (|6a|) . 

Parameters T and are related to (E^N) by the 
saddle-point conditions [1, Q 

g^E = 27r J dkuj^akbl , g^N = 2^ j dkauhl . (7) 

Given the saddle-point configuration one evalu- 

ates the suppression exponent. 



Fn{E) = g^ {2lmS[^s] - 2ET - NO) , 
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where the last two terms are due to non-trivial initial 
state, see Ref. 0. Note that (j)s{t,x) and Fn{E) depend 
on g, E, N only via combinations g'^N^ g^E^ see Eqs. ([7j). 

To summarize, the recipe for calculating the suppres- 
sion exponent of SA pair production in two-particle col- 
lisions is as follows. One starts at Sp > by finding two- 
parametric family of saddle-point configurations (j)s{t^x) 
which satisfy the classical field equation and boundary 
conditions (|6]). Then one computes the values of E^ N 
and Fn{E)^ Eqs. ([7j), (j8j), for each of these configurations. 
The result for the suppression exponent F{E) in Eq. ([1]) 
is recovered in the limits Sp ^ and g^N 0. 

We support the method by performing explicit calcu- 
lations in the model ([2]) with potential shown in the inset 
in Fig. [H 
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where W{x) = e~^^ {x + + x^), a = 0.4 and v is tuned 
to provide the required value of dp. We do not use the 
standard (j)^ potential because of the chaotic properties 
of (j)"^ kink-antikink dynamics [13, HI w hich lead to diffi- 
culties in the semiclassical analysis [ill, |l6| . 

We solve the semiclassical boundary value problem nu- 
merically [13 using methods of Refs. [3, Our start- 
ing point is bounce Q — Euclidean solution describing 
spontaneous decay of false vacuum at £^ = TV = 0. By 
using the classical field equation we continue the bounce 
to the Minkowski parts of the contour in Fig. ^jp. Then, 
changing (T, 0) (and thus {E^ N)) in small steps and solv- 
ing numerically the boundary value problem (j6j) , we con- 
struct the continuous family of saddle-point configura- 
tions (j)s{t^x) dit E < Ecb' Each configuration is repre- 
sented by a point in the left part of Fig. [TJ The points 
form the lines = const (filled points) and T = const 
(empty points). 

Solutions at £^ < Ecb^ Ecb ^ 2Ms have the form of 
distorted bounces, see Fig. [2t. Wave packets in the left 
part of the figure represent particles moving in the ini- 
tial state; after collision, the particles back-react on the 
Euclidean part of solution. Using the semiclassical solu- 
tions, we calculate the multiparticle exponent Fn{E) at 
E < Ecb- 

We evaluate the two-particle exponent F{E) by ex- 
trapolating Fn{E) [i9| to g'^N = with quadratic poly- 
nomials in g'^N^ cf. Eq. (jlj. The accuracy of extrapola- 
tion is 5%. 

The probability ([T]) of SA pair production must van- 
ish below the kinematic threshold E = 2Ms. Let us 
show that, indeed, F{E) +oo as ^ in the re- 
gion E < Ecb 2Ms. We recall that at small Sp the 
thin-wall approximation [5] is applicable and the two- 
particle exponent F{E) can be evaluated analytically 
modulo 0{Sp^) corrections [6], 
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We see that F{E) oc 1/Sp dit E < 2Ms; this property dis- 
appears for E > 2Ms where the thin-wall approximation 
breaks down. 

Our numerical results for 8p ■F{E) at g^N = (Fig. [3^, 
dashed lines) approach Eq. (p^Oj) (Fig. [3^, solid line) as 
Sp decreases and coincide with it after extrapolation to 
6p = (points). This gives support to our method. 

The fact that Sp ■ F{E) ^ as £; ^ 2Ms, Sp ^ {) 
hints that the properties of semiclassical solutions be- 
come qualitatively different E > Ecb- This is indeed 
the case; in fact, our numerical procedure does not pro- 
duce solutions at ^ > Ecb at all. By inspecting (j)s(t^x) 
with E ^ Ecb^ Fig. [2ll, we see the reason for that: this 
solution has long, almost static part where it is close to 
the critical bubble. The instability of the latter makes 
the numerical techniques inefficient. 
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Eq. (10) 
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FIG. 3. (a) Function Sp • F{E) at g^N = 0, ^ < 2Ms. (b) 
Suppression exponent Fn{E) ai E > 2Ms. 



We solve the instability problem by e-regularization 
method of Refs. [20|, EH which prescribes to add small 
imaginary term to the potential, V{(j)) Ve{(l)^x) = 
V{(j)) — ief{(j)^x)^ where f{(j)^x) > at real ^, 
/((/)+, x) = 0. For e > the critical bubble has com- 
plex energy and cannot be approached by (j)s{t,x): the 
energy of the latter is real due to Eq. (|6b|) . Thus, the 
semiclassical solutions are no longer unstable at e > 
and we are able to cover the region E > Ecb with new 
solutions, see Fig. [H Using Eq. (|8]) and taking the limit 
e ^ 0, we compute the multiparticle suppression expo- 
nent Fn{E). The typical saddle-point configuration at 
E > Ecb is shown in Fig. [2^. It still describes formation 
of the critical bubble [2l|; the energy excess E — Ecb is 
radiated away in the middle part of the solution. 

It is worth noting that formation of classically unstable 
"states" is a manifestation of general tunneling mecha- 
nism proposed recently in multidimensional quantum me- 
chanics TgI, [20I, [22I, [l2 and quantum field theory Til, [111 . 

To check the semiclassical method, we find values of 
{E^ N) corresponding to classically allowed decay of false 
vacuum 



24j . On the one hand, these values are obtained 
from the semiclassical exponent, since classical transi- 
tions correspond to Fn{E) = 0. Figure [T] shows the 
line Fn{E) = obtained in this way. On the other 
hand, "classically allowed" values of (E, N) can be ob- 
tained by studying classical solutions [25[. To this end, 
we performed [24] Monte Carlo simulation over the sets 
of Cauchy data {(/)(0, x), ^(0, x)}. Starting from each set, 
we solved numerically the classical field equation and ob- 
tained solution (j)ci{t^ x). We selected solutions describing 
transitions between the false and true vacua and com- 
puted the values of {E^N)^ Eqs. ([5]), ([7]), for these solu- 
tions. The latter values are shown in Fig.[T]by dots which 
fill precisely the region bounded by our line Fn{E) = 0. 
We conclude that the semiclassical result for Fn{E) is 
reliable. 

The exponent Fn{E) is plotted in Fig. [3b at A/" = 2/g^ 
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and different values of 5p (dashed lines). The graphs are 
almost indistinguishable; thus, the limit 5p ^ {) exists. 
Extrapolating Fn{E) to Sp = 0, we obtain the final result 
for the suppression exponent of SA pair production in N- 
particle collisions (solid line in Fig.lSb). The two-particle 
exponent F{E) is recovered by extrapolating Fn{E) to 
g^N = (upper graph in Fig. [Sb). 

At high energies our suppression exponents are almost 
constant which is a t ypica l behavior for collision-induced 
suppressions [26|, [Hf, Hsj . We expect that this feature 
holds in other models. 

We end up by arguing that the method proposed in 
this Letter is applicable in multidimensional {D > 2) and 
gauge theories. The problem with soliton-antisoliton at- 
traction can be solved in the general case by introducing 
small constant force F, an analog of Sp^ which drags soli- 
ton and antisoliton apart. At |F| 7^ SA pair is separated 
from the perturbative sector by a potential barrier; tun- 
neling through this barrier is a Schwinger process of spon- 
taneous soliton-antisoliton pair creation in the external 
field [23. For example, in the case of t'Hooft-Polyakov 
monopoles the external force is introduced by adding 
uniform magnetic field. Other steps of the method — 
RST conjecture (j4]) and semiclassical equations (j6j), ([7]), 
(jH]) — are explicitly general, cf. Ref. [H. In particu- 
lar, suppression exponent of Schwinger process isjro- 
portional to 1/|F| [27|; this enhancement vanishes [28] at 
E = 2Ms. The final result for the semiclassical expo- 
nent of collision-induced SA pair production should be 
obtained by taking the limit |F| ^ at > 2Ms. 

We thank V.A Rubakov, S.M. Sibiryakov, 
F.L. Bezrukov and G.I. Rubtsov for discussions. This 
work was supported in part by grants NS-5525.2010.2; 
MK-7748.2010.2, grant of the "Dynasty" foundation 
(D.L.); RFBR-ll-02-01528-a and Russian state con- 
tracts 02.740.11.0244, P520, P2598 (S.D.). Numerical 
calculations have been performed on the Computational 
cluster of the Theoretical division of INR RAS. 



* ldemidQv@ms 2.inr.ac.rul 
^ rievkov@ms2.i nF.ac.rul 
[1] A. K. Drukier and S. Nussinov, Phys. Rev. Lett. 49, 102 
(1982). 

[2] S. Dutta, D. A. Steer, and T. Vachaspati, Phys. Rev. 

Lett. 101, 121601 (2008); T. Romanczukiewicz and 

Y. Shnir, ibid. 105, 081601 (2010). 
[3] T. Banks, G. Farrar, M. Dine, D. Karabali, and B. Sakita, 

Nucl. Phys. B 347, 581 (1990); V. 1. Zakharov, Phys. 

Rev. Lett. 67, 3650 (1991); H. Goldberg, Phys. Lett. B 

257, 346 (1991). 
[4] Upon field rescaling (j) gcj), the coupling constant g 

enters the non-linear terms in the potential only, as it 

should. 



[5] M. B. Voloshin, 1. Y. Kobzarev, and L. B. Okun, Sov. J. 

Nucl. Phys. 20, 644 (1975); M. Stone, Phys. Lett. B 67, 

186 (1977); S. Coleman, Phys. Rev. D 15, 2929 (1977). 
[6] M. Voloshin and K. Selivanov, JETP Lett. 42, 422 

(1985); Sov. J. Nucl. Phys. 44, 868 (1986); V. A. 

Rubakov, D. T. Son, and P. G. Tinyakov, Phys. Lett. 

B 278, 279 (1992); V. G. Kiselev, Phys. Rev. D 45, 

2929 (1992). 

[7] A. N. Kuznetsov and P. G. Tinyakov, Phys. Rev. D 56, 
1156 (1997). 

[8] V. A. Rubakov, D. T. Son, and P. G. Tinyakov, Phys. 
Lett. B 287, 342 (1992); V. A. Rubakov and M. E. 
Shaposhnikov, Phys. Usp. 39, 461 (1996). 
[9] For alternative approach see S. V. lordanskii and 
L. P. Pitaevskii, Sov. Phys. JETP 49, 386 (1979); 
S. Yu. Khlebnikov, Phys. Lett. B 282, 459 (1992); D. Di- 
akonov, V. Petrov, Phys. Rev. D 50, 266 (1994). 
10] P. G. Tinyakov, Phys. Lett. B 284, 410 (1992); A. H. 

Mueller, Nucl. Phys. B 401, 93 (1993). 
11] G. F. Bonini, A. G. Cohen, C. Rebbi, and V. A. Rubakov, 

Phys. Rev. D 60, 076004 (1999). 
12] D. G. Levkov, A. G. Panin, and S. M. Sibiryakov, J. Phys. 

A: Math. Theor. 42, 205102 (2009). 
13] The saddle-point method is justified at small when 
the integrand exp(i5'[(/)]) in the path integral rapidly os- 
cillates, cf. Eq. dll). 
14] D. K. Campbell, J. F. Schonfeld, and C. A. Wingate, 

Physica D 9, 1 (1983). 
15] O. Bohigas, S. Tomsovic, and D. Ullmo, Phys. Rep. 223, 
43 (1993); A. Shudo and K. S. Ikeda, Phys. Rev. Lett. 
74, 682 (1995); 76, 4151 (1996). 
16] T. Onishi, A. Shudo, K. S. Ikeda, and K. Takahashi, 

Phys. Rev. E 64, 025201 (2001); 68, 056211 (2003). 
17] We discretize Eqs. (|6|), ([7|), (|8|) by introducing uniform 
lattice NtxNx ^ 5000 x 250 of extent Lt x = 15 x 7. 
Precision of discretization is kept smaller than 0.1%; it is 
controlled by changing the lattice size and keeping track 
of energy conservation. Linear evolution in Eq. (O holds 
with accuracy smaller than 1%. Extrapolation 6p ^ 
produces relative errors in F{E) of order 10~^. 
[18] F. Bezrukov, D. Levkov, C. Rebbi, V. Rubakov, and 
P. Tinyakov, Phys. Rev. D 68, 036005 (2003); Phys. 
Lett. B 574, 75 (2003). 



[19] One can argue M, \2$ that g^N 



as g^N 0, so 



we use only the first two terms in Eq. (|8|). 
[20] F. Bezrukov and D. Levkov, arXiv:quant- 
ph/0301022 (2003); JETP 98, 820 (2004); D. G. 
Levkov, A. G. Panin, and S. M. Sibiryakov, Phys. Rev. 
Lett. 99, 170407 (2007). 
[21] The critical bubble decays in the right part of Fig. [2^ due 

to regularization e > 0. 
[22] K. Takahashi and K. S. Ikeda, J. Phys. A: Math. Gen. 

36, 7953 (2003); FPL 71, 193 (2005). 
[23] D. G. Levkov and S. M. Sibiryakov, Phys. Rev. D 71, 

025001 (2005); JETP Lett. 81, 53 (2005). 
[24] S. V. Demidov and D. G. Levkov, JHEP 1106, 016 
(2011). 

[25] C. Rebbi and R. Singleton, Phys. Rev. D 54, 1020 (1996). 
[26] M. Maggiore and M. Shifman, Nucl. Phys. B 371, 177 
(1992). 

[27] N. S. Manton, Nucl. Phys. B 126, 525 (1977); 1. K. 

Affleck and N. S. Manton, 194, 38 (1982). 
[28] A. K. Monin, JHEP 0510, 109 (2005). 



